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Abstract. In this article we present the results from mid-infrared spectroscopy of a sample of 14 T Tauri stars with silicate 
emission. The qualitative analysis of the spectra reveals a correlation between the strength of the silicate feature and its shape 
similar to the one which was found recently for the more massive Herbig Ae/Be stars by van Boekel et al. (2003). The com- 
parison with theoretical spectra of amorphous olivine ([Mg,Fe]2Si04) with different grain sizes suggests that this correlation 
is indicating grain growth in the disks of T Tauri stars. Similar mechanisms of grain processing appear to be effective in both 
groups of young stars. 
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1. Introduction 

T Tauri stars are low mass young stellar objects with ages 
of 10^ - 10^ yr. One of the characteristics of so-called clas- 
sical T Tauri stars is a distinctive IR excess. The origin of 
the IR radiation is the accretion disk which plays a cru- 
cial role in star formation by transporting material to the 
central star while angular momentum is transferred out- 
wards. Different models were proposed to bring the observed 
facts into a general framework ( Lynden-Bell & Pringle ~974| 
lAdams, Lada, & ijhu l^tj7| [TFiang & Goldreich 1997^ Yet 
there remain many open questions concerning the details of 
configuration, chemistry and evolution of disks. A particu- 
lar feature of the infrared spectra of such disks is the sili- 
cate band in the region from 8 to 12//m. The band originates 
from the stretching mode of the Si-O bond of silicate min- 
erals like olivine ([Mg,Fe] 28104), forsterite (Mg2Si04), en- 
statite (MgSiOs) and silica (SiOa). The silicate feature ob- 
served in T Tauri stars appears in emission as well as in absorp- 
tion JCohen & Witteborn 19851 1. Silicate emission is assumed 
to emerge from a warm, optically thin disk layer (disk atmo- 
sphere) which is heated by the radiation of the central star 
( jChiang & Goldreich 1997] jNatta, Meyer, & BeckwithlOOO ). 
Honda et al. (2003) showed that crystalline silicates are present 
in T Tauri stars, indicating substantial grain processing. In this 
Paper we present studies on the silicate emission from a sam- 
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pie of 14 T Tauri stars based on mid-infrared spectroscopy per- 
formed during three observation campaigns in February, June 
and December 2002 on the ESO 3.6 m telescope at La Silla. 
We analyse the spectra in particular with respect to a possible 
correlation between shape and strength of the silicate feature. 
Such a correlation has been recently observed for Herbig Ae/Be 
stars ( van Boek el et al. 2003l l and has been interpreted as evi- 
dence of grain processing in the circumstellar disks of those 
stars. 

2. Observations and data reduction 

The observations of the T Tauri stars were performed with 
TIMMI2, the Thermal Infrai-ed Multi Mode Instrument 2 
( Reimann et al. 1998j installed at the 3.6 m telescope at ESO's 
observatory at La Silla, Chile. We observed each object in two 
modes: first we used the imaging mode to obtain photometry 
at 11.9yum (band width 1.2jum). In a number of cases in addi- 
tion we used the 8.9 fim filter (band width 0.8 /urn). Second we 
performed spectroscopy using the low-resolution grism mode 
(R=160) with a slit width of 1.2". In both modes the system 
was working with chopping and nodding with an amplitude 
of 10". In our data reduction, the measurements of the spec- 
trophotometric standards were used to determine the shape of 
the spectra, and the photometric data at 11.9//m were used 
to establish the absolute flux level. The spectrum of the at- 
mospheric extinction, which is needed to perform the first of 
these steps, was determined from observations of standard stars 
at different airmass. We regularly observed standard stars to 
monitor valuations in atmospheric ti-ansmission. The flux data 
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Table 1. List of observed sources with the siHcate feature in 
emission. 



Object 


Flux at 


Observed 


Multiple 


ISOPHOT 




n.9fim 


at Run " 


Source 


spectrum 




(Jy) 






available 


AK Sco 


2.1 


B 


* 


* 


AS 205 


8.4 


B 


* 




CR Cha 


0.8 


C 




* 


DRTau 


2.0 


C 




* 


GG Tau 


0.8 


C 


* 




GWOri 


6.3 


C 


* 




Glass I 


8.8 


A,B,C 


* 


* 


Haro 1-16 


0.9 


B 






HBC 639 


2.3 


A 






RULup 


2.4 


B 






SCrA 


5.1 


B 


* 


* 


SU Aur 


3.6 


C 


* 




TWHya 


0.6 


C 






WW Cha 


5.3 


B,C 







" A: 5. +6. Feb. 2002, B: 7.+8. Jun. 2002, C: 24.-27. Dec. 2002 
for AS 205 and S CrA we were able to obtain separate spectra of 
the components 



for the spectrophotometric standard stars, (based on models by 
ICohen 1998K were taken from ESO's TIMMI2 webpage. 

In total we observed 21 TTauri stars with the silicate feature 
in emission. We selected the 14 objects for which the silicate 
feature was measured with a signal-to-noise ratio better than 2 
(see Tab.^). Several of these are known to be binary or multiple 
systems. In two cases we were able to obtain separate spectra 
of the components. The remaining 7 sources (DO Tau, DQ Tau, 
FUOri, BBW76, VZCha, VWCha, Sz82) were degraded by 
noise to the extent that a meaningful spectral analysis was not 
possible. For some stars the observations were repeated to ver- 
ify the results. 

For our studies we assumed that the continuum outside of 
the silicate band is actually reached at the edges of our mea- 
sured spectral range (at 8.2/im and at 13.0yum). This as- 
sumption is supported by comparison of our spectra with 
ISOPHOT spectra (Abraham et. al 2003 in prep, see also 
|Natta , Meyer , & Beckwith 2 000 1. Figure ^ illustrates this for 
two examples. Here, the continuum was estimated by fitting a 
second order polynomial to all data points shortwards of 8 /im 
and beyond 12.8 /^m. Note that the fitted curve is matching the 
ends of the TIMMI2 spectra and has an almost linear charac- 
teristic over the range of the silicate feature. Since we were not 
able to find such additional data of good quality for all objects 
of our sample, we estimated the continuum level within the sil- 
icate feature by a linear interpolation between the end points of 
the TIMMI2 spectra at 8.2yum and at 13.0jum, where the values 
for these end points were obtained by averaging over the data 
over an interval of 0.2 yum. The continuum level determined this 
way agreed with the second order fit, where available, to typi- 
cally + 10 % . For all sources we consistently used the derived 
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Fig. 1. Comparison of TIMMI2 and ISOPHOT data for two ob- 
jects. The combined data were used to estimate the continuum 
by fitting a second order polynomial to the data points left of 8 
jum and longwards of 12.8 fim. 



linearly varying continuum level to determine the continuum 
normalized' spectra; these are shown in Fig.|2 



3. Analysis of the silicate feature 

The spectra in Fig. |2l show a variety in strength and shape. In 
some spectra (AKSco, Haro 1-16, CRCha) a strong silicate 
emission with a peak near 9.8/im is visible. In other cases 
(RuLup, DRTau, HBC 639), the silicate emission is weaker 
and the profile looks more like a plateau. It is known that 
the shape of the band is strongly affected by the chemical 
composition and in particular by the size of the dust grains 
(Henning et al. 1995, Bouw man et al. 2001 ). A coiTelation be- 
tween strength and shape of the silicate emission could be ex- 
pected if there are systematic changes in particle size within 
the group of objects. 

To estimate the strength of the silicate emission, we use the 
peak of the continuum normalised spectrum. For the charac- 
terisation of the band shape we measured the continuum sub- 
tracted flux at two wavelengths, 11. 3 and 9.8 fim. The flux ratio 
for these wavelengths can be interpreted as a indicator for the 
grain size composition following Bouwman et al. (2001). They 
point out that (see Fig.|3j 

- the mid-infrared absorption coefficient of small amor- 
phous olivine dust grains (with a size of 0. 1 /vm) has 
a triangular shape with the maximum at 9.8 jum. This 



to preserve the shape of the emission feature we define: 



^'norm(^) = 1 + [^totalC^) - F, 
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Fig. 2. Continuum normalized spectra of our sample ordered by the strength of the silicate feature. The shape of the feature is 
showing a correlation with the strength. Stronger features have a triangular shape with a pronounced peak near 9.8 fim while 
weaker features are more plateau-like. The gap from 9.15 to 9.65 /jm in most of the spectra is caused by a broken channel of the 
TIMMI2 detector. 
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Fig. 3. Theoretical spectra of olivine dust with grain sizes of 
0.1 yum and 2.0//m (from lBouwman et al. 20011 . The points at 
9.8 yum and 11.3 fim are visualizing the dependence of the flux 
ratio at these wavelengths on the grain size. 



spectrum is compatible with observations of the sili- 
cate emission of unprocessed dust of the ISM (see also 
|Bowey, Adamson, & Whi ttet 1998|l. 
- the mid-infrared absorption coefficient of larger amorphous 
olivine dust grains shows a plateau-like shape in the range 
from 9.5 to 12jum. 



Assuming a dust temperature of around 300 K, the blackbody 
thermal emission will not vary very much over our wavelength 
band, and the silicate emission profiles would have a shape 
quite close to the shape of the absorption coefficients. Since 
the absorption cross section per unit mass of large grains is in 
general smaller than that of small grains of similar shape, the 
emission of large grains is in general expected to be weaker 
than that of small grains. 

Our observations suggest that most observed mid-infrared 
spectra of young low or intermediate mass stars can be char- 
acterised simply by emission corresponding to the two grain 
sizes covered in Fig.|3 We note that the ratio of absorption co- 
efficients at the wavelengths of 11.3 fj.m and 9.8 jum as func- 
tion of grain size varies from 0.43 (for 0. 1 fj.m grains) to about 

I. (2.0jum grains). Therefore, the ratio of the (continuum sub- 
tracted) flux at 11.3jum over that at 9.8 yum should be a mea- 
sure of the grain size. In Fig.|4]we plotted this ratio against 
the strength of the silicate feature. The distribution of the data 
points shows a clear correlation. Strong silicate features have 

II. 3/9.8 flux ratios down to 0.53 (for CR Cha), while weak fea- 
tures exhibit higher ratios up to about 1 .0 (in case of the S CrA 
NW-component even 1.42). The remaining points are dis- 
tributed over the range between these two extremes. The error 
bars quoted in Fig.|4]take into account the signal to noise of the 
observations, and the uncertainty in the continuum estimate. A 
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0.4 [ ^ ^ ^ ^ - 

1 1.5 2 2.5 3 3.5 

feature strength [max (F^^,^)] 

Fig. 4. Correlation between the strength of the siHcate feature 
and the ratio of the (continuum subtracted) fluxes at 11.3/im 
over that at 9.8 fxm. The strength is defined by the maximum 
of the siHcate feature over the continuum. 

weighted linear fit to the data points with the relation y - a-bx, 
where x is the strength of the feature and y the flux ratio, 
gives the coeflicients: a - 1.11 +0.11 and b - -0.16 + 0.04. 
A very similar correlation between strength and shape of the 
silicate emission was found by van Boek el et al. 20031 for iso- 
lated Herbig Ae/Be stars, with coefficients of a = 1.48 + 0.09, 
b = -0.28 + 0.04. 

Emission by Polycyclic Aromatic Hydrocarbons (PAHs) at 
11.3 fim does not appear to be a problem for our data set. There 
was no evidence for PAH emission in our observed spectra 
nor in the available ISOPHOT-S spectra. This absence of PAH 
emission is not unexpected, since one can assume that the UV 
radiation of T Tauri stars normally is too weak to cause a sig- 
nificant excitation of PAH l |Natta & Kruegel 1995) . 

4. Discussion 

The large diff'erence in geometric cross section per unit mass 
between large and small dust grains implies that if the sili- 
cate emission band is dominated by small grains, a substantial 
amount of large grains could be present without being noticed. 
On the other hand, the spectra with plateau-like emission are 
in geneal weak, likely due to a much lower number of small 
grains present in the upper disk regions. In some spectra (S CrA 
both components, RULup) the flux at 11.3/im is even higher 
than the flux at 9.8/im, which results to a ratio greater than 
1. This cannot be explained by the assumption of amorphous 
oUvine dust with the absorption spectra shown above. It may be 
that crystalline forsterite (Mg2Si04) is responsible for this ef- 
fect, since its absorption spectrum shows a peak at 1 1 .3 yum (see 
[Bou wman et al. 2001 Hon da et al. 2003> . Possible indications 
for crystalline forsterite are also seen in some other objects with 
stronger features (SU Aur, GG Tau). 

We also note that the silicate emission emerges from warm, 
optically thin material, assumed here to be located in the up- 
per layer of a flared accretion disk, which is irradiated by the 
central star. The dust in these upper layers of different T Tauri 



disks appears to vary continously between two extreme types: 
a) small amorphous dust grains or a mixture of small and large 
grains, with a spectrum very similar to the dust of the ISM 
and b) large amorphous dust grains with a possible addition 
of crystalline silicate. If we think of these diff'erences in terms 
of evolution, then the reason for the difference between these 
two dust types could be either the removal of the small grains 
from the upper disk layers or the coagulation of small grains 
to larger ones. A removal of small dust grains by pure settling 
to the mid-plane of the disk is unlikely to explain the observed 
kind of differentiation since the settling time for larger grains is 
shorter than for the smaller ones. A further mechanism which 
could be responsible for the vanishing of small grains is by the 
radiation pressure of the central star A star of solar luminosity 
would be quite effective in removing tenth micron sized parti- 
cles. But in this case, a permanent diffusion of particles from 
the mid-plane on shorter time scales should avoid a significant 
decrease of the amount of small particles in the upper layer 
(ITakeuchi & Lin 2003J . The effect of coagulation is a possible 
explanation, if constructive collisions to larger grains dominate 
over destructive collisions which produce again small grains. 
This should be the case in TTauri disks where the SED can be 
well fit by gas-rich disks in hydrostatic equilibrium. 
An interesting detail of our studies is the availability of sep- 
arate spectra from the components of two binary T Tauri sys- 
tems (AS 205 and S CrA). We observed significant differences 
in strength and also differences in the shape of the silicate band 
of the components, with the secondary having the - relatively - 
stronger emission. Since the components presumably have the 
same age, the difference must have other reasons than the evo- 
lutionary state of the stars. 

Acknowledgements. This work is partly supported by the Hungarian 
Research Fund (OTKA No. T037508) and by the Bolyai Fellowship 
(P.Abraham). S.Y. Melnikov thanks the German Academic Exchange 
Service (DAAD) for the DAAD fellowship grant. 

References 

Adams, F C, Lada, C. J., & Shu, F H. 1987, ApJ, 312, 788 

van den Ancker, M. E., Wesselius, P. R., Tielens, A. G. G. M., van 

Dishoeck, E. F, & Spinoglio, L. 1999, A&A, 348, 877 
van Boekel, R., Waters, L. B. F. M., Dominik, C, Bouwman, J., de 

Koter, A., Dullemond, C. P, & Paresce, F 2003, A&A, 400, L21 
Bouwman, J., Meeus, G., de Koter, A., Hony, S., Dominik, C, & 

Waters, L. B. F M. 2001, A&A, 375, 950 
Bowey, J. E., Adamson, A. J., & Whittet, D. C. B. 1998, MNRAS, 

298, 131 

Chiang, E. I. & Goldreich, P 1997, ApJ, 490, 368 
Cohen, M. & Witteborn, F C. 1985, ApJ, 294, 345 
Cohen, M. 1998, AJ, 115, 2092 

Henning, T., Begemann, B., Mutschke, H., & Dorschner, J. 1995, 
A&AS, 112, 143 

Honda, M., Kataza, H., Okamoto, Y. K., Miyata, T., Yamashita, T., 

Sako, S., Takubo, S., & Onaka, T. 2003, ApJ , 585, L59 
Lynden-Bell, D. & Pringle, J. E. 1974, MNRAS, 168, 603 
Natta, A. & Kruegel, E. 1995, A&A , 302, 849 
Natta, A., Meyer, M. R., & Beckwith, S. V. W. 2000, ApJ, 534, 838 
Reimann, H., Weinert, U., & Wagner, S. 1998, Proc.SPIE, 3354, 865 
Takeuchi, T. & Lin, D. N. C. 2003, ApJ, 593, 524 



